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Abstract 

A method to treat the ground state correlations beyond the RPA is presented. A 
set of nonlinear equations taking into account effects of the ground state correla- 
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tions on the pairing and phonon-phonon coupling is derived. The influence of such 
correlations on properties of the vibrational states in spherical nuclei is studied. 
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1 Introduction. 



The Random Phase Approximation (RPA) is the basic method to describe 
the nuclear vibrational motion. It is well known that the RPA violates the 
Pauli principle and many attempts have been performed to improve it ([1- 
15]). There is a rather complete list of references on that subject in paper [10]. 
Renormalized RPA equations that include corrections for the ground state 
correlations (GSC), have been applied not only to study properties of the low- 
lying isoscalar vibrations in spherical nuclei [11,12], but also to investigate 
the charge-exchange modes in nuclei [16,17] and the giant resonances in metal 
clusters [18]. To describe excited states all above mentioned works use the one 
phonon states that can be created acting by the relevant collective operator 
on the ground state of the system, which is treated as the vacuum state. From 
another point of view it is well known that due to the anharmonicity of vibra- 
tions there is a coupling between one-phonon and more complex states [19,20]. 
Taking into account such a coupling it is possible to describe particularities 
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of the low-lying states and damping of the giant resonances [20]. Up to now 
such coupling was considered for the RPA phonons only [20]. 

In the present paper we use phonons of the extended RPA (ERPA) [12] as a 
basis on which the quasi-particle phonon model (QPM) [20] equations are gen- 
eralized so as to account for the GSC in the description of nuclear vibrational 
states constructed by one- and two-phonon configurations. Besides the GSC, 
we take into account the Pauli principle corrections arising in the two-phonon 
terms due to the fermion structure of the phonon operators. As an example 
we study the effect of the GSC on the energies, transition probabilities and 
transition densities of the low- lying vibrational states in 68 Zn and compare 
present results with the results within other approaches. 

2. Basic formulae 

We employ the QPM - Hamiltonian including an average nuclear field de- 
scribed as the Woods-Saxon potential, pairing interactions, the isoscalar and 
isovector particle-hole (p-h) residual forces in separable form with the Bohr- 
Mottelson radial dependence [19]: 

H = E{£(^ " K)a] m a jm - -G^ : (P t P ) r : - 

r jm 

\ £ (4 A) +^ A) ):(M] M (r)M v ( ( xr)):} (1) 

A/«T=±1 

We sum over the proton(p) and neutron(n) indexes and the notation {r = 
(n, p)} is used and a change r <-> — r means a change p <-> n. The single-particle 
states are specified by the quantum numbers (jm); Ej are the single-particle 
energies; A T is the chemical potential; and are the strengths in the 
p-p and in the p-h channel, respectively. The pair creation and the multipole 
operators entering the normal products in (1) are defined as follows: 

K = E(- 1 ) j_m « m (2) 

jm 

Mt (r) = -=L= £ (-I)^(jW - m' | \n)f^ m a fml (3) 

jj mm 

where fty are the single particle radial matrix elements of residual forces. 

In what follows we work in quasiparticle (qp) representation, defined by the 
canonical Bogoliubov transformation: 

4m = U J a jm + (-l) j_m W-m (4) 
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The Hamiltonian can be represented in terms of bifermion quasiparticle oper- 
ators (and their conjugate ones): 

B(jj;Xfi) = H-l/^'OWm I MHm^'-m' (5) 

mm 

^ + (jj';V) = ^(jmjV|A/i)«+ m «t m , (6) 

mm' 

The phonon creation operators are defined in the 2-qp space in a standard 
fashion: 

QU = \ E{^> Xfj) - {-l) x -»<p% A(jj'; A - (J,)} (7) 

33 

where the index A = 0, 1, 2, 3, ... denotes multipolarity and /x is its z-projection 
in the laboratory system. The following relation can be proved using the exact 
commutators of the fermion operators: 

(0 | [Qa^QJ^,,] I 0) = \8 X ^ £(1 -q 03 ^4 - ¥#¥#](8) 

33 



where | 0) is the phonon vacuum, q-y = qj + qy and qj is the quasiparticle 
distribution in the ground state: qj = (2j + 1)~^(0 | B(jj; 00) | 0). 

The phonon and the pairing characteristics are determined by the following 
non-linear system of equations: 
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The formulae for the quasiparticle energies Sj = y A% + — A r ) 2 and for 
the coefficients Uj, Vj remain the same as in the usual BCS theory; the new 

(r) 

values for A T = £(1 - 2g i )(2j + l) Uj v A and A T , come from the eqs. (12) 

j 

and (13); e^> = Sj + e^, uty = UjV^ + VjUy. The pairing vibrations (A = 0) 
have been considered in [23]. 

The system of nonlinear equations (9)-(14) includes effects of the isoscalar and 
isovector forces and it is a generalization of equations derived in [1,11,12,22]. 
This system treats the GSC self-consistently and describes the coupling be- 
tween different vibrations, different phonon roots of a certain multipolarity 
and the pairing. The present scheme is called Extended RPA (ERPA). The 
factors (1 — Qjj'), distinguishing the new equations from the conventional BCS 
and RPA ones, take into account the blocking effect due to the Pauli princi- 
ple. If we put qjj' = in the r.h.s. of eq. (14), we get the expression for the 
quasiparticle distribution in the ground state in the RPA case [3,9]. 

The GSC affect not only the RPA, but they also should change the quasiparticle- 
phonon coupling. To take into account such effects we follow the basic ideas 
of the QPM. Hereafter we derive the generalized QPM equations which take 
into account the GSC beyond the RPA. As it was shown in our previous paper 
[23] the pairing vibrations give a negligible contribution to qj. On the other 
hand the two-phonon configurations including the pairing vibration phonons 
have an energy essentially higher than the configurations constructed from 
usual vibration phonons. That is why we do not take into account the cou- 
pling with the pairing vibrations in what follows. Using the completeness and 
orthogonality conditions for the phonon operators one can express bifermion 
operators A + and A by phonons: 

A + (jj';X^) + A(jj';\-ri = 

(! - Q S f) + <t>%)(Qti + R A ^ Qx-*) (is) 

i 

The initial Hamiltonian (1) can be rewritten in terms of quasiparticle and 
phonon operators in following form: 

H = h + h pp + h QQ + h QB (16) 

h + h pp = J2 £ j a tm a jm (17) 
jm 
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E ^=((-) A ^V, + Qx-JBrUf, A - /x) + fc.c.(19) 
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where t> ./ 
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UjUji — VjVji and 



yXi = yXi + yXi 



1 - (4 A) + k[ x) ) X? 



(20) 



f21^ 



One can prove that the solutions of the system of equations (9)-(14) obey the 
following equality: 



(Qx^i I H | QXfii) = ^xi 



(22) 



The term Hqb is responsible for the mixing of the configurations and, therefore, 
for the description of many characteristics of the excited states of even-even 
nuclei. In the simplest case the wave functions of those states could be written 
down as: 



^(A/i) = {E^)^+ E PxSW [Qt^Qt^JW (23) 

i X\iiX2t2 



with the normalization condition: 



(tt„(JM) | * V (JM)) =Yl R i( J ») + 

i 

2 E (P A A -(^)) 2 (1 + ^ J (A^,A 2 . 2 )) = 1 

Al«lA2i2 



(24) 



where 



K J (X 1 i 1 ,X 2 i2) = K J (X 1 i 1 ,X 2 i2 \ X 2 i2,X 1 i 1 ) 



and 
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K J (X 2 i 2 , Xi | Xi, \ 2 i 2 ) = (2A + 1)(2A 2 + 1) 



1 
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The mean value of H is 
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(*„(JM) I H I *„(JM)) =Yl R KJv)"Ji + 2 E (P£S(Jv)) 

i AiiiA2i2 

(u Xlil + cj A2i2 + Acj J (Aiii, A 2 i 2 ))(l + ^ J (Aiii, X 2 i 2 )) + 
2 £ ^(^^(^^(^(l + ^Aiii.Aaia)) 



where the matrix element coupling one- and two-phonon configurations 



(Qju\h QB \ QUQ 



• A2 «2 



f )^(J l /)(l+^(A 1 i ll A 2 j 2 )) 



^(A0 = (-l) Al+A2+A ^(A0 
^(AO^E^fi(^r) 



(r) 



U\ltl{Xi, r) = ^=y/(2X 1 + 1)(2A 2 + 1) £ (1 - q hh ) x 



313233 
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I x Xlil + X Xli 
Aa; J (Aii 1? X 2 i 2 ) = --. Y\[ \ . =K J (X 2 i 2 , A^ \ A^, \ 2 i 2 ) + 

4 ir y/ypny™ 

X ^ 2 + X ^ K J (X 2 i,X 1 i 1 I A^A^)] (29) 

y/y^y™ 

Calculating the mean value of H we used so called quasidiagonal approxima- 
tion for the quantities K J (X 1 i 1 , \ 2 i 2 | A 3 i 3 , A 4 i 4 ) because the diagonal terms 
dominate over nondiagonal ones (see [20,21]). 

Using the variational principle in the form: 

5 ( ^(A/i) \H\ ^(A/i)> - E v {{ * v {\n) | ^(A/x)) - 1) = (30) 

one obtains the following system of equations: 

{u Ji -E v )R i {Jv)+ £ P^(Ju)U^(Ju)(l + K J (X 1 z 1 ,X 2 t 2 ))=0(31) 

Ai«iA 2 i2 

2(u Xlh + u X2l2 + A W J (A A) \ 2 i 2 ) - E v )P££(Jv) +Y,Ri(Jv)U$ 1 2 (Ji) = 0(32) 



The energies of the states (23) are solutions of 

Ut&MU^iXi 1 ) (l + K J (X 1 i 1 ,X 2 i 2 ] 



F(E V ) = det 



(uxi - E v )8i. 



Aiii,A2i2 



^Aiii + ^A 2 i 2 + Aw J (Aiii, X 2 i 2 ) - E v 



= 0(33) 



The rank of the determinant (33) is determined by the number of the one- 
phonon configurations included in the first term of the wave function (23). 

The above derived equations have the same form as the basic QPM-equations 
[20,21] and we call them the extended QPM with the Pauli principle correc- 
tions (EQPMPP) in what follows. The GSC affect phonon energies u>\i, nor- 
malization constants y x 1 and renormalize the matrix elements of the quasiparticle- 
phonon interaction. If we put K J '(A^i, X 2 i 2 ) = we get the equations already 
derived in [22] where all the fourth-order terms in phonon amplitudes ip Xl , and 

(f) Xl ., were neglected (we call this approach the extended QPM (EQPM)). In 
the limit q--> — > one reproduces precisely all the expressions of the QPM 
with taking into account the Pauli principle corrections [20,21] (QPMPP). In 
the case when q--i = and K J (X±ii, X 2 i 2 ) = we have equations describing 
coupling of one- and two- RPA phonons without taking into account the Pauli 
principle (QPM approach in following discussions). 
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3. Results and discussion 



As an example we have performed numerical calculations for the 68 Zn. The 
Woods-Saxon potential parameters in use are from [23], slightly modified to 
better describe the ground state density. In so far as our single-particle spec- 
trum includes the bound and quasibound states we do not use any effective 
charge to calculate the electromagnetic transition probabilities. The pairing 
constants are fixed so as to reproduce the odd-even mass difference of 
neighbouring nuclei. The strength parameters for the cases based on the 
RPA and ERPA schemes are adjusted so that the B(EX) values calculated 
with the wave function (23) within both approaches are reasonably close to 
the experimental ones. This means that for the ERPA calculations are 
larger than for the RPA ones [12]. No changes of have been done for cal- 
culations without the two-phonon terms. We use the ratio / Kq^ = —1.2 
that enables one to reproduce excitation energies of the isovector giant reso- 
nances in spherical nuclei. It is worth to mention that in our previous papers 
[11,12,23] we took into account only the isoscalar interaction for the p-h chan- 
nel, but the inclusion of the isovector interaction does not affect the structure 
of low-lying states practically. 

Our studies in Zn isotopes [12,23] of the effect of coupling of vibrations with 
different multipolarities (A) (via ) show that in realistic calculations one 
can keep only A=2 and 3 mainly. Solving nonlinear equations (9)-(14) one 
can find the phonon amplitudes,energies and the quasiparticle distributions 
within ERPA. Making use of them as input values it is possible to define from 
equations (24),(31)-(33) energies and the structure of the states described by 
the wave function (23). 

Knowing the wave functions it is not difficult to calculate any matrix ele- 
ments and physical quantities. For example, the charge transition density is 
calculated by the formula: 



Pi j \r) = T,PUr){ 1 -(l - q jlj2 )u% 2 Y.RiWWl* + & 



v 



3W2 

(-) 



3132 



E V(2Ai + l)(2A 2 + l)P££(Ji/)X;(l - Qjsn) * 



Alii A2«2 




.73 



}^3s&+^n&} (34) 



A2*2 j,Aiii> 



The expression for the two-quasiparticle transition density pj^Jj) can be 



found in[24]. Our charge transition densities are folded with the formfactor 
of the proton charge distribution [25]. Using quantities (34) one calculates the 
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reduced transition probabilities from the ground to the excited state (JV)[26] 

oo 

B(EJ; 0+ -> (Jv)) = (2J + 1) | | r J+2 pi J) (^)^ I' • (35) 

o 

The results of our calculations for the quasiparticle distribution in the ground 
state of 68 Zn are shown in the table 1. It contains the values of qj obtained 
in the RPA and ERPA schemes. As one can see from table 1 the qj have 
large values for the subshells near the Fermi surface only and the ERPA gives 
stronger correlations in comparison with the RPA. A similar behaviour of q/s 
has been found for other Zn isotopes [12]. This is valid for our choice of the 
multipole constants. As it was mentioned above the multipole constants 
have been chosen to describe with a reasonable accuracy experimental B(E\)- 
values. The value of nf ] = 0.0259 MeV/fm for the non-linear problem is quite 
larger than the critical RPA constant k 2) = 0.0242 MeV/fm where the RPA 
solution becomes complex. (In the RPA = 0.0227 MeV/fm). The 

octupole constant in use are equal to = 0.0235 MeV/fm for the RPA and 
k = 0.0250 MeV/fm for the ERPA calculations respectively. The last value is 
smaller than the critical RPA constant for the octupole vibrations in contrast 
to the quadrupole ones. All the calculations based on the RPA phonons (RPA, 
QPM, QPMPP) have been performed with the RPA constants as well all the 
calculations based on the ERPA phonons (ERPA, EQPM, EQPMPP) have 
been done with the ERPA set of constants. It is worth to mention that as 
in the case of metallic clusters ([18]) our ERPA calculation gives also weaker 
correlations compared to the RPA ones, if the same set of multipole constants 
is used in both approximations, but in this case neither energies nor transition 
probabilities can be reproduced within the ERPA or its modifications. 

To study the influence of the GSC on the quasiparticle-phonon coupling we 
calculated the structure of the low-lying states in 68 Zn with the wave function 
(23). Experimental data [27,28] and results of our calculations for the 2^ 2 and 
3f states within different approaches are shown in the table 2. 

One can see from table 2 that the RPA and ERPA overestimate the energies 
and fail to reproduce the transition probabilities for the 2% and 3j~ states. 
Taking into account coupling of the one- and two-phonon components im- 
proves essentially the description of all states under consideration. Besides 
the transitions to the ground state, one can reproduce the B(E2)-value for 
the E2-transition between the first and the second 2 + states. The inclusion 
of the Pauli principle corrections in two-phonon terms changes to worse the 
description of the second 2 + state mainly. The EQPMPP describes energies 
and transition probabilities better than the QPMPP. It is worth to point out 
that the B(E2)- values for the E2-transition between the first and the second 
2 + states depend essentially on the two-phonon components of the wave func- 
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tion of the 2^ state and the last ones can be affected by the three-phonon 
terms, which are out of the present consideration. One can conclude that the 
most consistent approach from theoretical point of view (EQPMPP) where 
the Pauli principle is taken into account in both one- and two-phonon terms, 
gives a rather good description of experimental data in general. 

The transition probabilities are the integral characteristics of the vibrational 
states and they are less sensitive to the details of the nuclear wave func- 
tions than the differential ones. As it was pointed out in our previous papers 
[11,12,23], GSC affect essentially the charge transition densities. Being the 
spatial overlap between the ground state wave function and the excited state 
wave function, the charge transition density provides a good test for nuclear 
models. The surface nature of the low-lying collective states predicted by cal- 
culations performed within the Hartree-Fock (HF) approach with effective 
forces [29] and the finite Fermi systems theory [30] has been demonstrated in 
the experiments on inelastic electron scattering from magic nuclei [31]. Exper- 
imental and theoretical (based on the random phase approximation (RPA)) 
studies of the charge transition densities [24,32] of the low-lying states in some 
spherical nuclei are in reasonable agreement, but the theory gives too large 
fluctuations of the transition densities in the interior region. All theoretical 
calculations in RPA, as in HF, give the same behaviour in the nuclear inte- 
rior, which indicates a systematic problem of a more fundamental nature (a 
detailed discussion can be found in refs. [33,34]). 

To test the developed approach we calculated the charge transition densities 
using eq.(34) in 68 Zn within different approaches. The figures 1 and 2 show 
the transition charge densities (p[f\r)) from the ground to the first 2 + state 
in 68 Zn. The experimental data [28] are presented as a shadowed area. Fig.l 
shows results of calculations based on the RPA phonons while calculations 
based on the ERPA phonons are presented in fig. 2. 

The RPA reproduces the behaviour of the charge transition densities qualita- 
tively, but it overestimates the interior part of the p\ (r). As one can see from 
fig.l, the inclusion of the two-phonon terms (see the QPM case) reduces the 
bump in the interior part of pf\r) by 17%. This is due to the reduction of the 
contribution of the first term in eq. (34) (i?j < 1). In spite of the dominance 
of the first term, the second one gives an additional reduction of p too. The 
Pauli principle correction for the two-phonon terms (QPMPP - case) change 
results not more than by 3%. 

Taking into account GSC beyond the RPA results in a suppression of interior 
oscillations by 9% in comparison with the RPA case (see fig. 2). Such a deple- 
tion is related with the Pauli blocking effect for the proton two-quasiparticle 
configuration {2p 3 / 2 , 2p 3 / 2 }, which is mainly responsible for the interior bump 
in the charge transition densities in the Zn isotopes. According to our RPA 
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calculations the proton two-quasiparticle configuration {2p 3 / 2 , 2p 3 / 2 } gives a 
contribution about 40% into the norm of the first quadrupole phonon in 68 Zn. 
The inclusion of the GSC redistributes the strength of this configuration over 
many phonon roots and as a result the contribution into the first root be- 
comes lower. The GSC suppresses the contribution of the partial proton two- 
quasiparticle {2p3/2,2p 3 / 2 } transition density for which q2 P3/2 has the biggest 
value and as it was mentioned above plays an essential role in the structure of 
the interior part of the transition density for the 2^ states. The configuration 
{I/5/2, 1/5/2} gives some contribution to the interior part too and the same 
mechanism of suppression takes place for it. It should be noted that the am- 
plitudes of the oscillations for the configurations with low orbital momenta are 
bigger than for the ones with high orbital momenta. That is because the single 
particle wave functions with low orbital momenta are mainly localized in the 
interior part of nuclei. Besides the blocking effect there are changes in coeffi- 
cients Uj*j 2 because of the influence of the GSC on pairing and in the phonon 
amplitudes. All these effects suppress the interior bump. Taking into account 
the phonon coupling we get an additional lowering of the interior bump in 
the charge transition density and a reason of such lowering is the same as it 
was discussed above for the coupling of the RPA phonons. Finally the EQPM 
approach gives 30% reduction of the (r) in the interior nucleus region. It is 
worth to note that the Pauli principle corrections in two-phonon terms change 
the results slightly, but such corrections must be taken into account because 
they are often responsible for the weak electromagnetic transitions between 
excited states. One can see from eq.(35) that the outer part of the charge 
transition density is responsible mainly for the value of the reduced transition 
probability. Since all above discussed approaches have very close values for 
this part of the charge transition density, the calculated B(E2) are close too. 

We would like to emphasize that the depletion effect discussed above can 
not be reproduced by any renormalization of the force strength kP^ in the 
RPA if one wants to describe energies and reduced transition probabilities. An 
influence of the p-p channel of the multipole forces on the charge transition 
densities is under an investigation now. As concerning the charge transition 
density for the second 2 + state, it is impossible to treat it without taking into 
account the three-phonon terms, which are out of the present consideration 
due to its numerical complexity. 

The charge transition densities for the 3j~ states in all Zn isotopes have a clear 
surface nature and there are no strong oscillations in the interior region of 
the nucleus because of a destructive interference of the 2-qp partial transition 
densities constructed from the single-particle wave functions with different 
parity. 

4. Conclusion 
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A consistent treatment of the ground state correlations beyond the RPA in- 
cluding their influence on the pairing and phonon-phonon coupling is pre- 
sented. A new general system nonlinear equations for the quasiparticle phonon 
model is derived. It is shown that this system contains as a particular case 
all equations derived for the QPM early. The system is solved numerically for 
first time in a realistic case for 68 Zn to study the effect of the GSC on the exci- 
tation energies, transition probabilities and charge transition densities of the 
vibrational states. Taking into account the GSC results in better agreement 
with experimental data for the characteristics of the low-lying states. 

Up to now all realistic calculations taking into account GSC beyond the RPA 
have been performed with separable effective forces only (see [11,12,16,17]) 
because even for such simplified forces it is not very easy to solve the set 
of nonlinear equations in a large configuration space. The finite rank sepa- 
rable approximation for the Skyrme forces suggested recently [35] opens a 
new possibility to treat GSC beyond RPA basing on more realistic nuclear 
interactions. 
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Table 1 

Quasiparticle distribution in the ground state of 68 Zn 







q r 


-lj 






neutrons 


protons 


nlj 


RPA 


ERPA 


RPA 


ERPA 


2si/ 2 


0.0047 


0.0134 


0.0122 


0.0189 


I/7/2 


0.0163 


0.0222 


0.0320 


0.0451 


2^3/2 


0.0540 


0.0651 


0.0862 


0.1053 


1/5/2 


0.0482 


0.0609 


0.0324 


0.0445 


2P1/2 


0.0899 


0.0994 


0.0359 


0.0461 


109/2 


0.0230 


0.0295 


0.0144 


0.0161 


2^5/2 


0.0072 


0.0087 


0.0042 


0.0047 
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0.030 - 




r (fm) 



Fig. 1. The transition charge density from the ground to the first 2 + state in 68 Zn. 
The dashed line corresponds to the RPA result; the dot-dashed line - QPM; the full 
line - QPMPP. Experimental data are presented by a shadowed area. 
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r (fm) 



Fig. 2. The transition charge density from the ground to the first 2 + state in 
68 Zn. The dotted line corresponds to the RPA result; the dashed line - ERPA; the 
dot-dashed line - EQPM; the full line - EQPMPP. Experimental data are presented 
by a shadowed area. 
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Table 2 

Energies and B(EA)-values for up-transitions to some first vibrational states of 68 Zn 



State 


2f 


2+ 


3f 




Energy 


B(E2) 


Energy 


B(E2) 


Energy 


B(E3) 






0+^2+ 




0+-2+ 
2+ -> 2+ 

z l ^ z 2 




0+ -► 3^ 


EXP. 


(MeV) 
1.077 


(e 2 fm 4 ) 
1266 


(MeV) 
1.883 


(e 2 fm 4 ) 
46 
287 


(MeV) 
2.751 


(e 2 fm 6 ) 
38400 


RPA 


1.360 


1290 


2.390 


75 


3.830 


48240 


QPM 


1.090 


1200 


2.060 


48 
231 


2.760 


38550 


QPMPP 


1.140 


1220 


2.180 


27 
351 


2.840 


39200 


ERPA 


1.330 


1250 


2.320 


106 


3.980 


43070 


EQPM 


1.080 


1170 


1.810 


47 
212. 


2.760 


34700 


EQPMPP 


1.080 


1270 


1.960 


38 
436 


2.750 


35660 
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